Topological color codes on Union Jack lattices: 
a stable implementation of the whole Clifford group 
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We study the error threshold of topological color codes on Union Jack lattices that allow for the full 
implementation of the whole Clifford group of quantum gates. After mapping the error-correction 
process onto a statistical mechanical random three-body Ising model on a Union Jack lattice, we 
compute its phase diagram in the temperature-disorder plane using Monte Carlo simulations. Sur- 
prisingly, topological color codes on Union Jack lattices have a similar error stability to color codes 
on triangular lattices, as well as to the Kitaev toric code. The enhanced computational capabilities 
of the topological color codes on Union Jack lattices with respect to triangular lattices and the toric 
code combined with the inherent robustness of this implementation show good prospects for future 
stable quantum computer implementations. 

PACS numbers: 03.67.Lx, 75.40.Mg, 03.67.Pp, 75.50.Lk 



I. INTRODUCTION 

Recently, the error threshold for topological color codes 
on triangular lattices has been computed using Monte 
Carlo methods [l[ , as well as using duality arguments ■ 
Understanding the error-correction properties of topolog- 
ical color codes is of paramount importance and serves 
as a benchmark for comparing to other quantum error- 
correction codes. Being one of the central figures of merit, 
error thresholds describe the tolerable value of noise be- 
low which a quantum code can perform quantum oper- 
ations without the impact of decoherence effects due to 
the environment. 

In particular, for topological stabilizer codes there are 
two main instances of quantum codes that are simple 
enough such that detailed studies of their computational 
capabilities can be carried out: the original toric code 
(TC) introduced by Kitaev Q and topological color codes 
(TCC) [J]. In both cases, quantum gates are imple- 
mented by external operations that act nontrivially on 
the degenerate ground-state manifold of a given quantum 
lattice Hamiltonian. The ground state harbors the topo- 
logical stabilizer code while excited states correspond to 
errors in the system. Within this scheme, topological 
quantum operations are performed based on ground-state 
properties [5[. Alternatively, it is possible to implement 
topological quantum gates based on quasiparticlc braid- 
ing operations @. 

In comparison to the toric code, TCCs can encode a 
larger number of qubits at a given surface of fixed topol- 
ogy and the variety of topological quantum gates that 
can be applied transversally is larger Q. In particular, 
it is the possibility to perform the whole Clifford group 
of quantum gates in a transversal and topological way 
that makes TCCs especially appealing in quantum infor- 
mation theory. The error threshold is one figure of merit 
for the performance of any topological stabilizer code. 



If this value is very small for a given code then all the 
advantages of the new code are meaningless in practi- 
cal applications because small amounts of external noise 
would spoil its stability. 

Improving the computational capabilities of a quan- 
tum information system generally comes at the price of 
decreased robustness against environmental noise. How- 
ever, the exact threshold value for the qubit error rate 
p in TCCs on triangular lattices has been computed nu- 
merically to be p c = 0.109(2) [l|, which agrees within 
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error bars with the value for the toric co de 
result — also confirmed by other studies @, [l2 T13f 
couraging: Topological color codes preserve a high error 
tolerance for quantum operations without performance 
trade-offs (i.e., it is possible to perform complex quan- 
tum operations such as quantum distillation, teleporta- 
tion, and dense coding while retaining the same stability 
against noise). 

To estimate the error tolerance of a TCC, the error- 
correction process is mapped onto a statistical model 
with random three-body interactions that correspond to 
the faulty bits. In analogy to the Kitaev model, the ran- 
dom three-body Ising model [l[ plays a similar role as the 
random bond Ising model in the Kitaev toric code [lij . 
The study of the three-body random Ising model further 
highlights the relationship between spin-glass physics and 
quantum information theory, and presents a new class of 
model systems exhibiting glassy behavior via three-body 
interactions (i.e., without spin- reversal symmetry). 

In this work we revisit the problem and estimate the 
error threshold for topological color codes on the Union 
Jack (UJ) lattice. The motivation to use the triangular 
lattice (which is dual to the hexagonal lattice) in Ref. [l[ 
was based on the triangular lattice being the simplest 
example of a family of lattices named colexes (color com- 
plexes) [HI and the simplicity of numerical simulations. 
However, there is a technical caveat with TCCs on the 
hexagonal lattice: they do not fully reproduce the Clif- 



ford group Q. In fact, one of the gate generators of the 
group needs the square-octagonal lattice for its topolog- 
ical implementation. We therefore use the fact that the 
UJ lattice is the dual lattice of the square-octagonal lat- 
tice. 

There is another fundamental reason to study TCCs 
on the UJ lattice. The presence of three-body random 
interactions poses new problems that were absent in the 
Kitaev code. For example, the type of lattice and its 
connectivity may play a role in some of the essential 
features of the mapped statistical model that is visible 
in the location of the multicritical point or the values 
of the critical exponents. To see this, there is another 
route to establish connections between quantum infor- 
mation and classical statistical mechanical models. It 
is based on the notion of classical simulability of topo- 
logical ground states [lf| [ISand classical models with 
the completeness property jl8l H^ |. The implication of 
these studies for TCCs is that certain properties of their 
ground state can be related to a three-body Ising model 
without randomness (p = 0) [^J. Following universality 
[20l l2l| , it is expected that the critical exponents depend 
on the lattice geometry and order parameter symmetry 
p2l I23I ]. Therefore, TCCs allow us to see a connection 
between critical exponents and computational capabili- 
ties of a quantum code, and not just the location of its 
critical point. This is a novel feature that is absent in the 
statistical mapping from topological codes and random 
models [l4| where only the location of the multicritical 
point plays a role in the features of the quantum infor- 
mation system. We conjecture that the same connection 
between lattice-dependent critical exponents and differ- 
ent quantum capabilities for TCCs holds in the presence 
of randomness p 7^ at the multicritical point . 

Our simulations show that the numerical value of the 
multicritical point p c for the three-body Ising model on 
random UJ lattices agrees within error bars with the 
value obtained for the triangular lattice [l| . In turn, this 
highlights the stability of TCCs to external noise. In 
addition, we compare the phase boundary between a fer- 
romagnetic ground state and a paramagnetic one for the 
three-body Ising model on both triangular (TR) and UJ 
lattices. Our results show that these are rather similar. 

Note that throughout this article we assume that ex- 
ternal error correction is carried out on the protected 
system represented by the topological color code. This 
is the standard notion of quantum error correction [2~H - 
[29[ | - Recently, experimental realizations of topological 
error correction have been implemented (30j . as well as 
experimental proposals for TCCs using Rydberg atoms 
[311 ] have been made. However, there are also more de- 
manding schemes in which the topological color code can 
be internally protected leading to the notion of a self- 
correcting quantum computer [32l HH . 

The paper is organized as follows: In Sec. lU we in- 
troduce topological color-code states on Union Jack lat- 
tices, followed by the mapping of the error-correction pro- 
cess of the TCC onto random three-body Ising models in 
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FIG. 1: Simulated Union Jack lattice: The sites at the edges 
of the squares (inset) have coordination 8, whereas the sites in 
the center of the square plaquettes have coordination 4. The 
Hamiltonian in Eq. ((2J is a sum over all possible triangles with 
a weight r multiplied by the product of the three spins S 1 , S 2 , 
and S 3 on the triangle's edges (inset). For the simulations, 
we use periodic boundary conditions (i.e., the lattice is placed 
on a torus). 

Sec. Mil In Sec. [IV] we describe our Monte Carlo simu- 
lations and how to locate the multicritical point p c that 
corresponds to the error threshold of the TCCs. Results 
from simulations are shown in Sec. [Vj followed by con- 
cluding remarks. 

II. COLOR CODES ON THE UNION JACK 
LATTICE 

A TCC can be obtained from any two-dimensional 
(2D) lattice in which all plaquettes are triangles and ver- 
tices are three-colorable such that no link connects ver- 
tices of the same color. It is also possible to work in the 
dual lattice (called a 2-colex, see Refs. [HE]) but here 
we prefer to work in the triangular lattice to have a more 
direct mapping, as done in Ref. [l|. The lattice is em- 
bedded in a compact surface of arbitrary topology. Since 
information is encoded in topological degrees of freedom, 
the code is nontrivial only when the topology of the sur- 
face is nontrivial (i.e., the genus g of the surface has to be 
g > 1). Note that data can also be encoded in a planar 
surface with holes and appropriate boundary conditions. 
This makes topolo gica l proposals more amenable to ex- 
perimental setups [3J, |35[. In this work we embed the 
system in a Union Jack lattice (see Fig. [TJ to be able to 
transversally implement all Clifford gates 

The construction of the code C is done in the following 
way: Consider a physical system with a qubit at each 
lattice triangle and introduce the following vertex opera- 
tors that generate the stabilizer group of C [28| . For each 
vertex v there are two types of operators that correspond 
to Pauli operators of X or Z type: 

X v := (g) X A , Z v := (g) Z A . (1) 

A:d£A A:tigA 

that is, a vertex operator acts on all nearby triangles 
which may be 4 or 8 (sec Fig. [1]). Vertex operators pair- 
wise commute and square to identity, so that they gen- 
erate an Abelian group S called the stabilizer. The code 
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C is defined as the subspace stabilized by iS and thus 
contains the states with X v = Z v = 1 Vw. 

What makes the UJ lattice special is that the vertex 
operators — the generators of the stabilizer — have support 
on a number of qubits that is a multiple of four. This 
fact, together with the properties of general color codes, 
gives rise to an important feature, namely that Clifford 
gates leave the stabilizer invariant. After the introduc- 
tion of suitable boundaries on the lattice, this allows for 
the transversal implementation of the Clifford group on 
color codes 



III. ERROR-CORRECTION MAPPING 

When an encoded state in a color code C is subject to 
errors, the first step for their possible correction is the 
measurement of vertex operators. The resulting collec- 
tion of ±1 eigenvalues is called the error syndrome and 
gives information about which errors occurred. Indeed, 
the errors need not be guessed exactly but only up to a 
stabilizer element. 

Color codes have a structure with stabilizer generators 
that are either products of X or Z Pauli operators, but 
not both. This allows us to treat bit-flip and phase errors 
separately: A-type (iT-type) errors produce violations of 
Z-type (X-type) vertex operators. The correction of each 
type of error can be expressed in homological terms (see 
Ref. [H) such that the syndrome gives the boundary of 
the error. In order to do a successful correction, errors 
must be guessed only up to homology. In contrast to 
toric codes where one deals with the usual homology of 
paths on a surface, in TCCs two different types (colors) 
of paths are allowed. 

Consider a standard error model based on stochastic 
errors in which phase errors Z and bit-flip errors X are 
uncorrelated and occur with probability p at each qubit. 
In topological codes, when p is below a threshold p c , er- 
ror correction can be performed with perfect accuracy 
in the limit of infinite system size. Above the thresh- 
old, error correction becomes useless in the same limit. 
In Ref. [lH this error threshold was computed for toric 
codes by mapping the error-correction process to a ran- 
dom bond Ising model. The corresponding mapping for 
color codes was carried out in Ref. p| , which showed that 
a random three-body Ising model appears. Note that 
other general extensions of statistical mappings are also 
possible in the context of topological subsystem codes 

Let us summarize what happens when we apply the 
mapping to the Union Jack lattice. First, the spins of 
the statistical model correspond to the stabilizer gener- 
ators. Thus, we have a classical spin at each vertex of 
the lattice. As for the three-body interactions, they cor- 
respond to the physical qubits and hence there is such a 
term per triangle. Each of these terms carries a random 



TABLE I: Simulation parameters: L is the system size, A sa 
is the number of disorder samples, t cq = 2 b is the number 
of equilibration sweeps, Tmin [Tmax] is the lowest [highest] 
temperature, and At the number of temperatures used. 



p 






L 


N sa , 


b 


T min 


T max 


JVt 


0.00 




12, 


18 


50 


18 


2.200 


2.350 


31 


0.00 




24, 


30 


50 


19 


2.200 


2.350 


31 


0.00 






36 


50 


20 


2.200 


2.350 


31 


0.01 




12, 


18 


5 000 


18 


1.900 


2.400 


51 


0.01 




24, 


30 


5 000 


19 


1.900 


2.400 


51 


0.01 






36 


5 000 


20 


1.900 


2.400 


51 


0.02 




12. 


18 


5 000 


18 


1.900 


2.400 


51 


0.02 




21. 


30 


5 000 


19 


1.900 


2.400 


51 


0.02 






36 


5 000 


20 


1.900 


2.400 


51 


0.03 




12. 


18 


5 000 


18 


1.700 


2.200 


51 


0.03 




21. 


30 


5 000 


19 


1.700 


2.200 


51 


0.03 






36 


5 000 


20 


1.700 


2.200 


51 


0.04 




12. 


18 


5 000 


18 


1.700 


2.200 


51 


0.04 




21. 


30 


5 000 


19 


1.700 


2.200 


51 


0.04 






36 


5 000 


20 


1.700 


2.200 


51 


0.06 




12. 


18 


5 000 


18 


1.600 


2.100 


51 


0.06 




21. 


30 


5 000 


19 


1.600 


2.100 


51 


0.06 






36 


5 000 


20 


1.600 


2.100 


51 


0.08 




12. 


18 


5 000 


18 


1.400 


2.000 


61 


0.08 




24, 


30 


5 000 


19 


1.400 


2.000 


61 


0.08 






36 


5 000 


20 


1.400 


2.000 


61 


0.10 - 


0.11 


12. 


18 


10 000 


18 


0.750 


2.600 


38 


0.10 - 


0.11 


24, 


30 


10 000 


19 


0.750 


2.600 


38 


0.10 - 


0.11 




36 


10 000 


20 


0.750 


2.600 


38 


0.10 - 


0.11 




12 


10 000 


22 


0.750 


2.600 


38 


0.12 




12. 


18 


5 000 


18 


0.750 


2.600 


38 


0.12 




24, 


30 


5 000 


19 


0.750 


2.600 


38 


0.12 






36 


5 000 


20 


0.750 


2.600 


38 



sign yielding the Hamiltonian 

H = -JJ2^S 1 A S 2 A S 3 A , (2) 

A 

where the coupling J is positive and the sum is over the 
product of the spins at the vertices of all triangles on the 
UJ lattice and ta = ±1 [37j |. The (classical) Ising spins 
S A can have the values ±1. The sign of the coupling 
constants r are independent random quenched variables, 
that are negative with probability p. 

At low temperature T and small disorder p, the sys- 
tem is ferromagnetic. Above a critical line p cr ; t (T), the 
model undergoes a phase transition to a paramagnetic 
state. The error threshold p c is given by the crossing 
point between p cr n (T) and the Nishimori line [III [3^ 

exp(-2 J/T) = (3) 



IV. NUMERICAL DETAILS 

When studying phase transitions of statistical systems, 
it is most favorable to study quantities whose finite-size 
scaling form is dimensionless (i.e., there is not a system- 
size dependent prcfactor in front of the scaling function). 
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Typical quantities are the Binder ratio [40( as well as the 
dimensionless two-point finite-size correlation length di- 
vided by the system size [4l[ . The latter has been shown 
to be a more robust measure of transition temperatures 
for systems without spin-reversal symmetry |l| because 
the Binder ratio becomes negative and steep at the tran- 
sition (i.e., pinpointing the critical temperature is diffi- 
cult). 

The transitions to a ferromagnetic phase are deter- 
mined by a finite-size scaling of the dimensionless two- 
point finite-size correlation length divided by the system 
size. We start by determining the wave- vector-dependent 
susceptibility 



i) 



(4) 



In Eq. ((3]) (• • • )t represents a thermal average and R; 
the spatial location of the spins. The correlation length 
is then given by 



1 



2sin(/c min /2) y [x(fcmin)]av 



1, 



(5) 



where fc m in = (2tt/L, 0) is the smallest nonzero wave vec- 
tor and [■•■]av represents an average over the different 
error configurations (disorder sampling). The finite-size 
correlation length divided by the system size has the fol- 
lowing finite-size scaling form: 



U/L-XiL^^T-Tc]) 



(6) 



where v is a critical exponent and T c represents the tran- 
sition temperature. Numerically, finite systems of linear 
size L are studied. In that case the function £ m /L is in- 
dependent of L whenever T = T c as then the argument of 
the function X is zero. In other words, when different sys- 
tem sizes arc studied, the data cross at one point, which 
corresponds to the transition, if present. This can be seen 
in Figl^l Because finite-size scaling corrections are small 
in this case, one can use the estimate of T c obtained as 
a very good approximation to the thermodynamic limit 
value. The critical exponent v for the correlation length 
can be determined by a full scaling of the data [lU, as 
shown in Ref. [![. 

We have also computed the spin-glass finite-size corre- 
lation length [replace Si with SfSf in Eq. ©]. For all 
values of p studied, no sign of a finite-temperature spin- 
glass transition was found. This result coincides with the 
results found for the triangular lattice [l[ and is compat- 
ible with the standard belief [39[ that no glassy phase 
is expected in two-dimensional random models, although 
here this is extended to systems with three-body interac- 
tions. We believe that a spin- glass phase may be possible 
for spin models with random many-body interactions in 
three space dimensions. 

Because the complexity of the problem increases con- 
siderable when p > 0, we use parallel tempering Monte 



Carlo [43|, |44|. In addition to local spin flips, after each 
lattice sweep a global update that exchanges the temper- 
ature of two replicas (copies) of the system is proposed. 
This considerably speeds up the simulations (parameters 
are shown in Table [I]). Equilibration is tested by a loga- 
rithmic binning of the data. Once at least the last three 
bins agree within errors, we define the system to be equi- 
librated. Note that we use periodic boundary conditions 
to reduce finite-size effects. To prevent a mismatch with 
the vertex coloring rules of the lattice and boundary con- 
ditions, we use system sizes L that are a multiple of 6. 



V. RESULTS 

Figure [2] shows the temperature-dependent dimension- 
less finite-size correlation length for different values of 
p. The top- left panel shows data for p = 0. The vertical 
dashed line is the analytically known transition tempera- 
ture of the two-dimensional Ising model (i.e., T c ~ 2.2692 
[2(il|). It is remarkable that the three-body Ising model 
on the Union Jack lattice has the same transition tem- 
perature, albeit being in a different universality class. Its 
critical exponents are fully determined by knowing two 
of them (e.g., v = 3/4 and a = 1/2 [H). For p w 0.108 
critical behavior sets in, and for p = 0.110 no sign of a 
transition is visible. We thus estimate p c = 0.109(2), in 
agreement with results for the triangular (TR) lattice [l[ 
and the toric code [8T4Tl|. 

The full p -T c phase diagram for the U J lattice is shown 
in Fig. [3J the solid (blue) line is the Nishimori line. To 
check the differences between the phase diagram for the 
UJ and TR lattices, we plot the relative deviation be- 
tween the different estimates of the critical temperature, 
(T C UJ - T c TR )/T c TR , as a function of p (Fig. gj. For most 
values of p the fluctuations are statistical and not larger 
than 1%. For p ~ p c , fluctuations are larger and the de- 
viations arc of the order of ~ 7%. Therefore, we believe 
that the phase boundaries for both models are very close. 



VI. CONCLUSIONS 

In this work we have computed the p-T phase dia- 
gram for the random three-body Ising model on a two- 
dimensional Union Jack lattice (Figs. [T] and [3]). The origi- 
nal motivation for this numerical study with Monte Carlo 
methods is to compute the multicritical point p c in this 
phase diagram: the crossing point between the critical 
line separating ordered/disordered phases and the Nishi- 
mori line. The crossing point p c corresponds to the er- 
ror threshold for topological color codes defined on the 
square-octagonal lattice (the dual of which is the Union 
Jack lattice) — a value that decides whether a topological 
stabilizer code is good enough for performing quantum 
error correction in practical applications. Our numerical 
result of p c = 0.109(2) is in agreement with the corre- 
sponding value for the TCCs on triangular lattices and 
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FIG. 2: (Color online) Finite-size correlation length £. m /L as a function of temperature T for different values of p. The data 
for p = cross at the critical temperature of the 2D Ising model (dashed line top left panel). For p < p c = 0.109 there is 
signature of a transition (data for different L cross) whereas for p = 0.110 > p c the transition vanishes. Note that for p > 0.108, 
corrections are large thus making the determination of the transition temperature difficult. For the case of p — 0.110 we only 
show the two largest sizes for clarity. 
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p c = 0.109(2) 



the toric code, within error bars. This result shows that 
TCCs on the square-octagonal lattice, which allow for the 
implementation of the complete Clifford group of quan- 
tum gates [H , are similarly stable as the toric code. The 
fact that the triangular lattice and the UJ lattice share 
similar values of p c that agree within error bars is by no 
means obvious since both models are known to be in dif- 
ferent universality classes. Note that the Monte Carlo 
method combined with finite-size scaling used to deter- 
mine the error threshold does not rely on approximations. 
The computed values of the phase boundary can merely 
be affected by finite system sizes, although corrections to 
scaling are small. Furthermore, effects of realistic error 
models have yet to be tested. 



0.02 0.04 



0.06 
P 



0.08 0.1 0.12 



FIG. 3: (Color online) p-T c phase diagram for the random 
three-body Ising model on the Union Jack lattice. For p > 
p c ~ 0.109 the ferromagnetic order is lost. The dotted line is 
a guide to the eye, the black circle represents the analytically 
known transition temperature of the 2D Ising model. The 
blue (solid) line represents the Nishimori line. In the regime 
marked by a dashed line the exact determination of T c (p) is 
difficult. 
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P 

FIG. 4: (Color online) Relative difference between T c for the 
Union Jack (UJ) and triangular lattices (TR) [l|] as a func- 
tion of p. In most cases the deviations are smaller than 1%. 
Around p c « 0.109 strong fluctuations appear since it is hard 
to estimate the transition temperatures. The data for the tri- 
angular lattice have been adapted from Ref. p} and expanded. 
Note that for p — > the deviations are smallest, suggesting 
that the slopes of the phase lines for p = are very close. 
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